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This article presents a set of curves that can be used to evaluate in detail 
over a wide range of operating conditions the performance of systems 
using envelope detectors to search for CW signals in narrow-band Gaussian 
noise. The charts discussed relate the probability, P D , of detecting such 
signals to (i) the signal-power /noise-power ratio, S/N, (ii) the proportion of 
false detections, P FA , and (Hi) the number, M, of independent samples 
of the envelope of the combined signal and noise that are averaged in making 
one attempt at detection. The curves and scales were calculated and drawn 
entirely by computer. The computation program was designed so that 
nearly linear curves coidd be produced, thereby increasing the accuracy 
and ease of interpolation. 

I. INTRODUCTION 

The curves* given in this article relate statistical properties of the 
output of an envelope detector to the signal-power/noise-power ratio, 
S/N, of a CW signal in narrow-band Gaussian noise at the input. 

The signal will be represented here by a finite section of a sine wave 
which may be divided further into segments, for convenience in proc- 
essing, so that most of the signal energy will be contained in a rela- 
tively narrow band, comparable in width to the reciprocal of the seg- 
ment length. A filter will be used to select the part of the spectrum in 
which the signal may be found. The output of this filter will be ap- 
plied to an envelope detector, and the amplitude of the output of the 
detector will be measured. Only if the magnitude measured at this 
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point exceeds a chosen threshold value will a signal be presumed 
present. 

For other detection tests two or more samples of the output of the 
detector, taken far enough apart in time that the noise values are in- 
dependent, will be averaged. Only if the value of such an average 
exceeds a chosen threshold will a signal be presumed present. The 
curves given here for the averages of multiple samples depend on the 
assumption that a detected signal remains steady while the samples 
are taken, and that the noise background is weakly stationary. 

This form of detection criterion has been discussed by several au- 
thors, 1 - 2 - 3 - 4 who give major emphasis to the use of a square-law de- 
tector because it is easier to analyze. A linear detector, however, is 
much more practical and is extensively used, so results that specifically 
apply to it are desirable. 

The curves available in the literature are difficult to use for checking 
system performance at more than a few scattered points. Those given 
here, on the other hand, are believed to be suitable for defining system 
performance in detail over a considerable range of threshold values 
and, furthermore, are easy to use. 



II. DISCUSSION OF THE MATHEMATICAL MODEL 

We assume that the spectrum of noise in the vicinity of the signal 
is smooth enough so that the output of the filter can be considered a 
narrowband Gaussian process with the possible addition of a sinusoidal 
signal. S. 0. Rice 5 - a derived an expression for the probability density 
of the envelope of such a waveform. Rice's expression, paraphrased, is 

P(p) dp = 2h 2 p exp [-h 2 (r 2 + p 2 )]I (2h 2 rp) dp (1) 

which gives the probability that the envelope lies within an interval, 
dp, of p. P(p) is thus the probability density. I (z) is the modified 
Bessel function of zero order defined by 

/.« = ± If- (2) 

Note that in (1) the term h 2 r 2 , which represents the signal/noise power 
ratio, and 1/ \/2h, which represents the rms noise level, can be identified. 
In this article S/N always means the ratio of the signal power to the 
noise power accompanying it in the specified narrow band. 
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When there is no signal (1) becomes 

P(p) dp = 2A a pexp (-h 2 P 2 ) dp (3) 

which represents a Rayleigh distribution having 

2a 2 = 1/h 2 (4) 

for the second moment about zero. 

The distribution of the envelope of the filter output is thus given by 
(1) when a signal is present, and by (3) when there is no signal. From 
these expressions distribution curves can be drawn for cases where 
only noise is present, and where a known signal occurs with the noise. 
To generate the performance charts given here, one has to be able to 
calculate the area lying above some threshold for distribution curves 
corresponding to no signal, and for those corresponding to all the values 
of S/N of interest. 

Only in the case of noise alone, for single-sample detection, is it 
possible to get an explicit formula for the area under portions of such 
a curve. Moments of the distribution given in (1) can be calculated 
with the help of Ref. 7. We give here the nth moment about zero: 

Vn = P P(p) P n dp 

= 2 exp (-&V) P /t 2 p n+1 exp (-h 2 p 2 )I (2h 2 rp) dp (5) 

= ( V2,)"r(^),/.-,(-|;l;-/ 1 v). 
In this h'-r- is the signal/noise power and 

2a 2 - 1/h 2 

as stated earlier. 

T(z) is the Gamma function and ,Fj (a;b;z) is the Confluent Hy- 
pergeometric function defined by 

(6; 

a(a + l)(a + 2) z 3 
+ b(b + 1)(6 + 2)'3-2-l "*""•• 
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When the signal/noise power becomes zero we have 

,„. = (V2.)"r(^). (7) 

This gives the nth moment about zero for the distribution due to noise 
alone. We therefore have 



■w-*.,F.(-|;i;-»v) 



(8) 



which gives the nth zero moment of the distribution when S/N = h 2 r 2 . 
We can thus derive a type-A Gram-Charlier 8 (G-C) series approxima- 
tion to the distribution curve for any S/N, and, since the series can be 
integrated, the area under any desired part of the curve can be found. 
Knowing the moments of the primary distribution, we can calculate 
the moments of a secondary distribution formed by averaging M sam- 
ples of the primary distribution. We thus form a G-C series which can 
be integrated to get the area under any desired part of the secondary 
distribution curve. 

Charts are included that cover cases in which, when only noise is 
present, the probability is as low as 10"* that the threshold will be ex- 
ceeded, even with the use of single samples in detection attempts. 

Producing a G-C series that would give results within 1 percent of 
the true value when the threshold is great enough to give such a low 
probability required evaluating many more moments than is usual in 
G-C series. In fact, 31 moments were calculated, which resulted in a 
series that oscillates about the true value of the distribution as the 
threshold is varied. A method that sometimes reduces these oscillations 
is to use special groupings of the terms, forming a so-called Edgeworth 
series. 8 It occurred to the author, however, that a better fit might be 
obtained by averaging the approximations given by several G-C series 
carried out to different numbers of terms since their oscillations as 
the threshold varies are not in phase. This involves very little more 
work than deriving the series that has the most terms, since it merely 
requires that the terms be weighted to represent the number of times 
they occur in forming the average. Fig. 1 shows the results obtained by 
forming G-C series that use up to 31 moments and averaging over the 
highest 23 and 13 terms, respectively. Both show less oscillation over 
restricted — though appreciable— ranges of threshold than does the 
Edgeworth series that uses up to 30 moments. The curve showing 
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THRESHOLD — *■ 

Fig. 1 — Comparison of methods for evaluating G-C series. 

the logarithm of Q, the function approximated by the series, allows 
the useful range to be assessed. It makes possible the calibration of the 
horizontal scale, which is the same for all four curves, in terms of 
the probability that the threshold will be exceeded. The logarithm to 
the base ten of this probability is given by the right scale of Fig. 1. 

III. DESCRIPTION OF THE CHARTS 

We will now describe briefly how the charts (Figs. 2 through 15) 
were produced and point out their main features. Some of these are 
believed novel and show the advantage of being able to use a digital 
computer to generate and draw charts of this type. 

For each value of S/N at which a curve was desired, the detection 
threshold, T, was expressed in two ways: (i) measured with respect 
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to the mean of the distribution for noise only, and standardized by 
dividing by the square root of the variance, and (ii) measured with 
respect to the mean of the distribution appropriate to the S/N, and 
standardized with respect to this distribution. 

The form in ({) is appropriate for calculating the probability of 
false alarm (i.e., the probability that noise alone will cause the thresh- 
old to be exceeded) , and the form in (ii) is appropriate for calculating 
the probability of detection (the probability that a signal at the chosen 
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Fig. 2 — ROC curves. 



LINEAR DETECTOR 



761 



0.300- • • ■ 
0.400- • • • 
0.500- • ■ • 


2 ir- 

S/N F 


D 


EP 

M 


ENDENT SA 
4.0 TO 12.0 


UPLES AVERAGE 
IN ldB INCREM 



ENTS 


\: . 




Jff.. ■ -■ • 

* .. * .... ' 


0.700- • ■ • 

> jS 

0.800 ««• ■ 




JT K 


0.900 •• - - 






* 












0.950 ^- • 
0.975 • 'S* 

0.995-/"— 









-,-5 5 2 , n - 6 5 



3 4 5 

STANDARDIZED THRESHOLD (UNIT IS «T M ) 

Fig. 3 — ROC curves. 



S/N will cause the threshold to be exceeded). The relationship between 
these two forms is linear, so that if the threshold were varied a straight 
line would give the relationship between a scale for measuring the 
variation based on (i) and another based on (ii). If these scales were 
calibrated in terms of probability of false alarm, P FA , and probability 
of detection, P D , useful charts could be obtained that use these prob- 
ability parameters as horizontal and vertical axes, respectively. 
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But, for any P PA scale, the corresponding Pj> scale depends on S/N. 
Fortunately, the P D scale does not change rapidly with S/N for the 
range we are interested in, and, to allow several S/N values to be in- 
cluded on each chart, a constant normal probability scale was chosen 
for P D . The penalty incurred by doing this is the slight curvature of 
some of the S/N lines plotted in Figs. 2 through 15. 

It can be seen that the horizontal scale, Pfa , changes as the num- 
ber of samples averaged changes. This is because the probability dis- 
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Fig. 4 — ROC curves. 
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Fig. 5 — ROC curves. 



tribution for the average of a number of samples depends on the num- 
ber averaged as well as on the basic distribution. Since an explicit 
formula for the cumulative probability is known in the case of the 
basic (Rayleigh) distribution, the Rayleigh distribution was used to 
get the Pfa scale for the chart based on single-sample decisions. All 
the other scales were calculated from G-C series approximations to the 
appropriate distributions. 



764 



THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1907 



To make interpolation easier, a column of heavy dots was put near 
the right edge of each chart. These dots mark 0.2 dB increments in 
S/N between the lines. Since any two adjacent S/N lines are very 
nearly straight and parallel it is easy to interpolate to within 0.1 dB 
between the lines using a parallel ruler guided by the heavy dots. 

The horizontal and vertical grids were plotted as rows and columns 
of dots in order to carry the accuracy of the scales into all parts of 
the chart without also producing a confusing mesh of lines. 
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Fig. 6 — ROC curves. 
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Fig. 7 — ROC curves. 



Each decade of the horizontal, P FA , scale, reading from left to right, 
is given by 

10.0(1.0)4.0(0.5)2.0(0.2)1.0, 



where the increments are in parentheses. 
The vertical, Pp , scale, reading from top to bottom, is 
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0.05(0.01)0.10(0.02)0.90(0.01)0.950(0.005)0.980 
0.980 (0.0025) 0.990 (0.001 ) 0.998 (0.0005) 0.999, 

the increments again being in parentheses. 

On Fig. 16 are curves for two pairs of P D , Pfa values that show how 
the S/N for each pair varies with the logarithm to the base 2 of the 



O.050 



64 INDEPENDENT SAMPLES AVERAGED 
S/N FROM -4.0 TO 1.0 IN I dB INCREMENTS 




3 4 5 

STANDARDIZED THRESHOLD (UNIT IS (T M ) 

Fig. 8 — ROC curves. 
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Fig. 9 — ROC curves. 



number of samples averaged. Over any restricted range in which the 
number of samples is doubled, either curve could be approximated 
quite satisfactorily by a straight line. 

The slopes of these two curves are almost the same at the same S/N, 
although the P It , P PA values are quite different. This means that, 
within any small range of Pp , Pfa values, curves like those shown in 
Fig. 16 can be assumed parallel with little sacrifice in accuracy. Con- 
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sequently, the change in S/N corresponding to a small change in log 2 Af 
indicated by such a curve will hold quite accurately for a small range 
of P D , P FA values around the pair for which the curve is actually 
drawn. 

The S/N corresponding to any pair of P D , P F a values can therefore 
be calculated for a different number of samples averaged by assuming 
a linear relationship between S/N and the logarithm to the base 2 of 
the number of samples. 
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Fig. 10 — ROC curves. 
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Fig. 11 — ROC curves. 



IV. CONCLUSION 



4.1 Efficiency of Post-Detector Integration 

When more than one sample can be averaged to improve detection 
sensitivity the process is usually called post-detector integration. The 
curves of Fig. 16 show how efficient this process is at different values 
of S/N and P D , Pfa operating points. The gradients of these curves 
show the rate at which detectability is improved, measured in dB/ 
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double-the-nuniber-of-samples, ds, for two operating points. This ratio 
is convenient to use in interpolating between adjacent charts, since the 
number of samples doubles from one chart to the next. 

It can be seen that the gradients are almost the same where the 
S/N is the same. The rate of improvement is better than 2 dB/ds for 
S/N above 10 dB, presumably reaching a maximum of 3 dB/ds for 
high S/N. The rate apparently becomes constant at 1.5 dB/ds for S/N 
values less than —10 dB. 
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Fig. 12 — ROC curves. 
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Fig. 13 — ROC curves. 



4.2 Extending the Scope of the Charts 

The detector input is restricted to a relatively narrow band of fre- 
quencies to ensure that the envelope can be measured accurately after 
detection. By bandshifting high enough, this result can be obtained 
even for bands of quite appreciable width. The charts can therefore 
be used in situations where the significant signals may be quite short 
bursts of nearly sinusoidal waveform. 
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The same standardized threshold scale applies to all the charts. It 
is given in units of 



a.\t = a/s/M, 



(9) 



where a is the rms noise into the detector and M independent samples 
are averaged. 

Although the curves given here refer strictly to envelope detectors 
they can be applied when square-law-detectors are used if a small 
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Fig. 14 — ROC curves. 
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Fig. 15 — ROC curves. 



error (always less than 0.2 dB) in S/X is taken into account, and the 
.standardized threshold scale is changed. A calibration curve for this 
error is given in Ref. 1. Reasonable agreement with this curve was 
found by generating ROC curves for a square-law detector after suit- 
ably modifying the computer program. It was possible thus to get 
ROC curves for 1 sample and for averaging 128 or more independent 
samples. 
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Fig. 16 — Sensitivity related to number of samples averaged. 
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